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An air bubble trapped in water by an oscillating acous- 
tic field undergoes either spherical or nonspherical pulsations 
depending on the strength of the forcing pressure. Two dif- 
ferent instability mechanisms (the Rayleigh-Taylor instability 
and parametric instability) cause deviations from sphericity. 
Distinguishing these mechanisms allows explanation of many 
features of recent experiments on single bubble sonolumines- 
cence, and suggests methods for finding sonoluminescence in 
different parameter regimes. 



The discovery and elucidation of single bubble sono- 
luminescence (SL) [1-3] sparked a renewal of interest in 
the dynamics of a levitated bubble driven by an oscil- 
lating acoustic field. Detailed optical measurements [4,5] 
(cf. figure 2 of [4]) reveal different dynamic regimes, dis- 
tinguished by the bubble shape, which depend on the 
ambient bubble radius Ro and the strength of the forc- 
ing pressure amplitude P a . For very small forcing, the 
bubble is spherical throughout the oscillation period. At 
higher forcing the bubble develops non-spherical surface 
oscillations. Two different instability mechanisms are op- 
erating: (i) The Rayleigh-Taylor (RT) instability [6], oc- 
curring whenever gas is strongly accelerated into a liquid, 
(ii) Parametric instability, arising due to the accumula- 
tion of perturbations from sphericity over many oscilla- 
tion periods. 

However, experiments show several phenomona which 
are confusing in light of these two instabilities: (a) Near 
the onset of sonoluminescence for air bubbles in water, 
the nonspherical pulsations vanish [1,4,5] and the bub- 
ble regains spherical symmetry. This is at variance with 
the RT instability, which should be more potent at the 
stronger accelerations required for SL. (b) In highly vis- 
cous fluids, the parametric instability is diminished. Nev- 
ertheless, it has proven extremely difficult (if not impos- 
sible) to make bubbles in highly viscous fluids (greater 
than about ten times the viscosity of water) [1] undergo 
strong oscillations, even at very strong forcing. 

The purpose of this paper is to assist in resolving these 
puzzles by presenting a stability analysis [7,8] of the os- 
cillating bubble. Our main result is that we account 
for both above experimental phenomena (a) and (b) by 
taking into consideration the detailed interplay between 
the two instabilities (i) and (ii). All parameters corre- 
spond to the experimental ones [3,4,9] for an air bubble 



in water: the surface tension of the air- water interface is 
a = 0.073kg / s 2 , the viscosity of water is v = 10~ 2 cm 2 / s, 
its density is p w = \000kg jrr? , the speed of sound in wa- 
ter is c w = 1481m/s, the driving frequency of the forcing 
acoustic field P(t) = P a cos(uit) is ui/2tt = 26AkHz, the 
external pressure Po = latm, and finally j = 1.4 for the 
ratio of the specific heats. 

The radius R(t) of a driven spherical bubble obeys the 
Rayleigh-Plesset (RP) equation [10,9,8] 
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Lofstedt, Barber, and Putterman show that if p(R,i) is 
adiabatically slaved to the bubble radius via a van der 
Waals equation of state, i.e. 

m = p °jihw> (2) 

then the Rayleigh-Plesset equation (1) well describes the 
experimental R(t) [9]. Here, h = i?o/8.54 (for air) is 
the hard core van der Waals radius. The adiabatic RP 
equation applies whenever R < c gas , where c gas is the 
speed of sound in the gas [9]. 

We now focus on the stability of the radial solution 
R(t). Following previous authors [7,11,8], consider a 
small distortion of the spherical interface R(t) + a n (t)Y n 
where Y n is a spherical harmonic of degree n. The dy- 
namics for the distortion amplitude a n (t) is given by 

a n + B n (t)a n - A n (t)a n = 0. (3) 

The full expressions for A n (t) and B n (t) for a gas bubble 
in a viscous fluid were first derived by Prosperetti [8]; 
they are nonlocal in time, reflecting the interaction be- 
tween the bubble and the flow in the fluid, initiated by 
the bubble wall motion itself. Our calculations employ a 
local approximation of the full hydrodynamic equations 
(17), (23), and (25) of [8], assuming that there will be 
fluid flow only in a thin layer of thickness 6 around the 
bubble. Approximate solutions of the full hydrodynamic 
equations [12] give the penetration depth 6 as the mini- 
mum of the diffusion length scale y 'v /us and of R/(2n). 
Then 

a / 1 \ / fin*? 

A " (t) = (n - l) R-p-R3 
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- ((" - 1)(" + 2) + 2n(n + 2)(n - l)-£) (4) 

and 

B„(i) = 3| + ((n + 2)(2n + 1) - 2n(n + 2) 2 -£) ^, 

(5) 

where f3 n = (n — l)(n + l)(n + 2). Note that the second 
term in eq. (5) causes damping [13] of the shape oscilla- 
tion and is always positive. The penetration depth takes 
into account both radial and azimuthal diffusive fluxes in 
the velocity field. The azimuthal gradients ensure that 
the viscous contribution to B n (t) stabilize the bubble. 

Rayleigh-Taylor Instability: First we focus on the RT 
instability, which occurs near the minimum bubble radius 
where the acceleration R of the bubble wall is positive 
(i.e., the gas accelerates into the fluid). Nonspherical 
perturbations of the bubble shape grow during this time 
period. The amplification factor follows from a WKB 
type analysis on equation (3), taking a n (t) ~ exp (S n (t)). 
The average amplification can be estimated as 



S n « J S n (t)dt = J (-^ ±yJ?i + A n jdt (6) 

where the integral is evaluated over the time period dur- 
ing which S n (t) is positive. Appreciable growth occurs 
when molecular fluctuation (size lnm) can grow to the 
minimum size of the bubble. Figure 1 shows the phase 
diagram of this instability as a function of the ambient 
bubble radius Ro and the forcing pressure P a . 

Parametric Instability: After transient effects, R(t) 
and thus A n (t) and B n (t) are periodic with period T = 
27r/w [14]. Eq. (3) is then called Hill's equation. Para- 
metric instability corresponds to a net growth of a non- 
spherical perturbations each oscillation period, so that 
after many periods perturbations overwhelm the bub- 
ble. Formally, this occurs whenever the magnitude of 
the maximal eigenvalue of the Floquet transition matrix 
F n (T) is larger than one. F n (T) is defined by 

(;;5)--m(«i(o))- 

By numerically computing F n (T) and determining its 
eigenvalues, we mapped out the entire phase diagram of 
stability. We first calculate the stability diagram for the 
n = 2 mode at zero viscosity (fig. 2). In general, when 
P a and Ro are large, the bubble becomes more unstable, 
although the detailed structure of the phase diagram is 
quite complicated. Many features can be understood an- 
alytically [11] by examining the small forcing (P a <C jPo) 
limit, where eq. (3) can be reduced to Mathieu's equa- 
tion. Substituting b n (t) oc (R(t)fl 2 a n (t) one obtains 



b^ + 2£ M b' n + lu 2 m (1 + e M cos(t))b n = 0, (8) 

with uj 2 M = l3 n a/(p w Rlu> 2 ), e M = P a /{'(Po) and £ M = 
[(n + 2)(2n + l)-2n(n + 2) 2 <V R ]v / (R 2 uj). Primes denote 
derivation with respect to the dimensionless time t = tot. 
The well known Mathieu tongues [15] occur in fig. 2 at 
u> M = k/2, where k is an integer. Since u>m oc Rq 3 ^ 2 
the tongues are not equidistant but become more and 
more packed at small Ro- For n = 2 the first Mathieu 
instability (k = 1) occurs at Ro = 50//m. 

Finite viscosity £m oc v ^ stabilizes the surface 
dynamics. Though is small (^m ~ 10 -4 for Ro = 
10/um), the Mathieu tongues are stabilized [16]. The sta- 
bilization is stronger at small Ro, where u>m and are 
large. Figure 1 shows a superposition of the stability di- 
agrams for modes n = 2,3,4,5,6 corresponding to the 
viscosity of water v = 10 _2 cm 2 / s. Stability diagrams of 
this type were first considered by Eller and Crum [11], 
and later by Horsburgh and Holt [17]. These studies 
examine larger bubble sizes than shown in fig. 1; our 
calculations in the large bubble regime [12] give similar 
thresholds as found in [11,17]. A particularly interesting 
new feature of our fig. 1 is the presence of small islands 
of stability dispersed throughout the unstable domain. 

When the bubble oscillations are weak, there is satu- 
ration of the linear instability, leading to an oscillating 
nonspherical bubble. Such shapes have been observed 
[11,17]. In the strongly nonlinear regime, saturation 
is unlikely, since parametric instability causes a n (t) of 
roughly constant amplitude throughout a cycle, whereas 
the bubble size R(t) changes by orders of magnitude. 
The RT instability causes 10 4 growth of perturbations 
on timescales of less than 10 _9 s, which almost certainly 
leads to the destruction of the bubble. 

What are the consequences of these results for the SL 
experiment? A standard experimental protocol [18,1,4] is 
to slowly increase the driving pressure P a for a bubble of 
given ambient radius Ro- For low viscosity fluids (c.f. fig. 
1) the parametric instability sets in before RT. The para- 
metric instability acts over a long time scale (~ 10 cycles) 
which is comparable to the diffusive timescale. The bub- 
ble therefore has time to readjust its size Ro by enhancing 
diffusion through the nonspherical surface , thereby reen- 
tering a stable parameter region [19]. In this way the long 
(10 _3 s) time scale parametric instability "protects" the 
bubble from encountering a region of the phase diagram 
where it could be destroyed by the (10 _9 s) RT instability. 
Upon further increasing P a , the same mechanism works 
again: The bubble more or less tracks the parametric in- 
stability borderline in the phase diagram until finally the 
parameter regime is reached where the forcing is strong 
enough so that a shock can be emitted from the collapsing 
bubble wall. The shock is associated with the sonolumi- 
nescence [21] and experimentally starts at P a ps 1.15atm 
[4]. In the strong forcing regime there may be additional 
mechanisms associated with the sonoluminescence itself 
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that help stabilize the bubble [22]. 

An immediate consequence of above scenario is that 
the bubble dynamics should be hysteretic: On increasing 
the forcing pressure, the bubble undergoes nonspherical 
oscillations. At an even slightly higher forcing pressure, 
the bubble will shrink. If the forcing pressure is then 
slowly decreased, the oscillations will become spherical 
again, even beyond that pressure where nonsphericity 
initially appeared. This type of hysteretic behavior was 
noted in the experiments of Gaitan [1]. 

We next discuss how increasing the viscosity changes 
the above scenario. A large viscosity increases the para- 
metric instability thresholds: at a viscosity of ten times 
that of water, the entire parameter range shown in figure 
(1) is stable. On the other hand, the Rayleigh-Taylor in- 
stability depends only weakly on viscous effects, because 
the maximum acceleration of the bubble is essentially in- 
dependent of viscosity [9]. Our numerical calculations 
show the RT threshold changes only slightly. Thus in 
highly viscous fluids (as glycerol), the bubble encounters 
the RT instability at a smaller driving pressure than the 
parametric instability, rather than the other way round, 
as in water. The RT instability acts so quickly (10 4 am- 
plification within 10 _9 s) that the bubble has no time 
to reenter a stable part of parameter space before being 
completely destroyed. We hypothesize that this could 
explain why SL has not yet been observed in highly vis- 
cous fluids. However, note that this problem could be 
surmounted by preparing or choosing bubbles of small 
enough ambient radius that as the driving pressure is 
continuously increased, the bubble does not experience 
the Rayleigh Taylor instability before sonoluminescing. 
An alternative explanation for the absence SL in glycerol 
relies on the very different solubility of air in water and 
glycerol, respectively [20]; we cannot rule this out. 

We again point out that our calculations depend on the 
basic assumption that the pressure is adiabaticity slaved 
to the bubble radius. This only holds when the bubble 
wall velocity is subsonic. For an air bubble in the Ro 
range of interest, adiabatic calculations show that the 
bubble wall becomes supersonic at P a ~ 1.35atm. The 
experiments [4,23] show that sonoluminescence begins al- 
ready around P a ~ 1.15atm (indicating the onset of some 
shock [21]). We take this as a hint that a more complete 
treatment should take heat diffusion [24,9] into account, 
since decreasing j decreases the speed of sound in the 
gas. Work is in progress [12]. 

In conclusion, we have analyzed the stability of a cav- 
itating bubble to both Rayleigh Taylor and parametric 
instabilities. At the high forcing pressures necessary for 
SL, only a bubble with a very small ambient radius is 
stable to both the RT and parametric instability. For 
low viscosity fluids, the parametric instability "pushes" 
the bubble to this special region by forcing the ambient 
size of the bubble to decrease. For more viscous fluids the 
parametric instability does not exist in the appropriate 



parameter region, so the bubble is destroyed by the RT 
instability before sonoluminescing. At the heart of the 
argument are the very different timescales on which the 
RT (10 _9 s) and the parametric instability (10 _3 s) act. 
Finally, our scenario also accounts for hystersis. 

Many mysteries remain: Are there additional mecha- 
nisms for stabilizing the bubble wall which arise from the 
sonoluminescence itself? What is the nature of the ambi- 
ent radius dynamics at a fixed forcing pressure? What is 
the role of mass transport mechanisms beyond diffusion 
[20]? Why are bubbles with inert gases so much stabler 
than nitrogen bubbles? 
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FIG. 1. Phase diagram for finite viscosity v = 10~ 2 cm 2 /s. 
Ro is given in /tin and P a in atm. The black region in- 
dicates stability with respect to all modes n = 2,3, 4, 5, 6. 
The diagram is strongly dominated by the instability of the 
n = 2-mode and considering even higher modes (n > 6) does 
not change the figure. Also shown is the RT stability curve 
(solid line) beyond which the bubble is RT unstable. Param- 
eters are given in the text. 

FIG. 2. The phase diagram of steady state bubble motions 
for n = 2 and v = 0. All other parameters are as in fig. 
1. The black region indicates stability, and the white region 
instability. 
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